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Abstract 

Extending the approach developed by Araiijo and Stoeger [l(] and improved in Araiijo et al 
4], we have shown how to construct dust-filled A 7^ Friedniann-Lemaitre- Robertson- Walker 
(FLRW) cosmological models from FLRW cosmological data on our past light cone. Apart from 
being of interest in its own right - demonstrating how such data fully determines the models - it 
is also illustrated in the flat case how the more general spherically symmetric (SS) Einstein field 
equations can be integrated in observational coordinates with data fit to FLRW forms arrayed 
on our past light cone, thus showing how such data determines a FLRW universe - which is not 
a priori obvious. It is also shown how to integrate these exact SS equations, in cases where the 
data are not FLRW, and the space-time is not known to be flat. It is essential for both flat and 
non-flat cases to have data giving the maximum of the observer area (angular-diameter) distance, 
and the redshift Zmax at which that occurs. This enables the determination of the vacuum-energy 
density //a, which would otherwise remain undetermined. 

PACS numbers: 98.80.-k, 98.80.Es, 98.80.Jk, 95.36. +x 



I. INTRODUCTION 

The recent WMAP's results [see Spergel et al [3|] and references therein] strongly sup 
port the inflationary scenario and are consistent with a nearly flat Friedmann-Lemaitre- 
Robertson- Walker Universe with a cosmological constant A and dust {Vt ^ 1, with Q^n ~ 0.27 
and ^A ~ 0.73), and with an almost scale invariant spectrum for the primordial perturba- 
tions. It is obviously very important to test and conflrm this result. One observationally 
sensitive theoretical approach to doing so begins by establishing a more general framework 
than FLRW - say a general perturbed spherically symmetric space-time - and then using 
the data itself to determine the more speciflc model. Can we obtain perturbed FLRW by 
doing this? 

In two papers Araiijo and Stoeger 1] and Araiijo, Roveda and Stoeger 2| demonstrated in 
detail how to solve exactly the spherically symmetric (SS) Einstein fleld equations for dust in 
observational coordinates without assuming FLRW and with cosmological data representing 
galaxy redshifts, observer area distances and galaxy number counts as functions of redshift. 
These data are given, not on a space-like surface of constant time, but rather on our past 
light cone C~{pq), which is centered at our observational position pq "here and now" on our 
world line C. These results demonstrate how cosmologically relevant astronomical data can 
be used to determine the space-time structure of the universe - the cosmological model which 
best flts it. This has been the aim of a series of papers going back to the Physics Reports 
paper by Ellis et al. [4]. The motivation and history of this " observational cosmology (OC) 
program" is summarized in Araiijo and Stoeger [l|. All these papers assumed that A = 0. 
In this paper we demonstrate how this program may be carried out when A 7^ 0. As a 
simple, and very relevant example, we take a flat SS universe, and suppose that the redshift, 
observer area-distance, and number-count data can be flt to FLRW functional forms (these 
are very special forms the data must take, if the universe is FLRW). Then we show how 
such data determines a bona fide FLRW universe - which is not a priori obvious. We then 
go on to indicate how the solution can be obtained in the more general non-flat case, with 
no constrains on the functional form of the data functions. In doing this for A 7^ 0, we 
also need another necessary piece of data, the maximum of the observer area distance, and 
the redshift at which it occurs. Without these extra observables, we do not have enough 
independent data to determine the model - in particular to determine the extra parameter 



A. 

The primary aim of OC program is to strengthen the connections between astronomical 
observations and cosmological theory. We do this by allowing observational data to deter- 
mine the geometry of spacetime as much as possible, without relying on a priori assumptions 
more than is necessary or justified. Basically, we want to find out not only how far our ob- 
servable universe is from being isotropic and spatially homogeneous (that, is describable by 
a„ FLRW cosmological model) on various length scales, but also to g,ve a dynamic account 
of those deviations (Stoeger et a/ [5[). 

By using observational coordinates, we can thus formulate Einstein's equations in a way 
which reflects both the geodesic flow of the cosmological fluid and the null geometry of 
C (Po)) along which practically all of our information about the distant reaches of our 
universe comes to us - in photons. In this formulation the field equations split naturally 
into two sets, as can be easily seen: a set of equations which can be solved on C^{po), that 
is on our past light cone, specified by w = wq, where w is the observational time coordinate; 
and a second set which evolves these solutions off C^{po) to other light cones into the past 
or into the future. Solution to the first set is directly determined from the data, and those 
solutions constitute the "initial conditions" for the solution of the second set. 

There are many reasons for investigating FRLW in observational coordinates from this 
more general starting point. It is clear from the Cosmic Background Radiation (CBR) 
anisotropics measured by the Cosmic Backgroung Explorer (COBE) and by WMAP that 
the universe is not exactly described by the FLRW models (this follows from the analysis 
by Sachs and Wolfe [a], and see Stoeger, Araiijo and Gebbie 71] for an analysis related to 
the viewpoint of this paper). But on the largest scales its deviations from FLRW are small. 
So on those scales the universe can be described by an almost FLRW model. Therefore, our 
first step towards a strictly observationally-based approach to this realistic model involves a 
complete understanding of the inner workings of the integration procedure in observational 
coordinates for FLRW data. 

In this paper, for completeness, we review some aspects of the problem of determining 
the solution of the exact spherically symmetric Einstein equations for dust in observational 
coordinates and then integrate the field equations with FLRW data to obtain the FLRW 
{k = and A 7^ 0) solution explicitly. We refer the reader to Ellis et al [4|, Kristian 
and Sachs [Sj, Araiijo and Stoeger [l| and references therein for a complete account of the 



philosophy and the foundations of the OC approach leading to the integration of Einstein 
field equations in observational coordinates. 

In the next section we define observational coordinates, write the general spherically 
symmetric metric using them and present the very important central conditions for the 
metric variables. Section IIIII summarizes the basic observational parameters we shall be 
using and presents several key relationships among the metric variables. Section HVl presents 
the full set of field equations for the spherically symmetric case, with dust and with A 7^ 0. 
Section |V] shows the integration procedure for the flat case, with FLRW data. In Section 
Wl\ we present the integration for the general, non-fiat case, and in section IVIII we briefly 
discuss our conclusions. 

II. THE SPHERICALLY SYMMETRIC METRIC IN OBSERVATIONAL COOR- 
DINATES 
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We are using observational coordinates (which were first suggested by Temple 
described by Ellis el al [4] the observational coordinates x* = {w, y, 9, cp} are centered on 
the observer's world line C and defined in the following way: 

(i) w is constant on each past light cone along C, with u°'daW > along C, where u"" is the 
4- velocity of matter {u°'Ua = — 1). In other words, each w = constant specifies a past light 
cone along C. Our past light cone is designated as w = Wq. 

(ii) y is the null radial coordinate. It measures distance down the null geodesies - with afiine 
parameter u - generating each past light cone centered on C. y = on C and dy/du > 
on each null cone - so that y increases as one moves down a past light cone away from C. 

(iii) 6 and are the latitude and longitude of observation, respectively - spherical coordi- 
nates based on a parallelly propagated orthonormal tetrad along C, and defined away from 
C by k"'da9 = k"'da4> = 0, where /c" is the past-directed wave vector of photons {k°'ka = 0). 

There are certain freedoms in the specification of these observational coordinates. In w there 
is the remaining freedom to specify w along our world line C. Once specified there it is fixed 



for all other world lines. There is considerable freedom in the choice oi y - there are a large 
variety of possible choices for this coordinate - the afiine parameter, z, the area distance 
C{w,y) itself. We normally choose y to be comoving with the fluid, that is u"'day = 0. Once 
we have made this choice, there is still a little bit of freedom left in y, which we shall use 
below. The remaining freedom in the 6 and coordinates is a rigid rotation at one point 
onC. 

In observational coordinates the Spherically Symmetric metric takes the general form: 

ds'^ = -A{w, yfdw"^ + 2A{w, y)B{w, y)dwdy + C{w, yfdn'^, (1) 

where we assume that y is comoving with the fluid, so that the fluid 4- velocity is u"" = A'^6^. 
The remaining coordinate freedom which preserves the observational form of the metric 
is a scaling of w and of y: 

f dw dy\ 

W^W = w{w), y^y = y{y) I — ^0^ — 1. (2) 

The first, as we mentioned above, corresponds to a freedom to choose w as any time 
parameter we wish along C, along our world line at y = 0. This is usually effected by 
choosing A{w,0). The second corresponds to the freedom to choose y as any null distance 
parameter on an initial light cone - typically our light cone a.t w = wq. Then that choice 
is effectively dragged onto other light cones by the fluid flow - y is comoving with the fluid 
4-velocity, as we have already indicated. We shall use this freedom to choose y by setting: 

Aiwo,y) = B{wo,y). (3) 

We should carefully note here that setting A{w, y) = B{w, y) off our past light cone w = wq 
is too restrictive. 

In general, these freedoms in w and y imply the metric scalings: 

A-.i=^A, B^B=^Ib. (4) 

dw dy 

It is important to specify the central conditions for the metric variables A{w, y), B{w, y) 
and C{w,y) in equation ([T]) - that is, their proper behavior as they approach y = 0. These 
are: 

as y ^0 : A{w, y) -^ A{w, 0) ^ 0, 



C{w,y)^B{w,Q)y = Q, (5) 

III. THE BASIC OBSERVATIONAL QUANTITIES 

The basic observable quantities on C are the following: 

(i) Redshift. The redshift z at time wq on C for a coninioving source a null radial distance 

y down C^(po) is given by 

/l(wo,0) 

This is just the observed redshift, which is directly determined by source spectra, once they 
are corrected for the Doppler shift due to local motions. 

(ii) Observer Area Distance. The observer area distance, often written as ro, measured 
at time Wq on C for a source at a null radial distance y is simply given by 

ro = C{wo,y), (7) 

provided the central condition ([S]), determining the relation between C{w,y) and B{w,y) 
for small values of y, holds. This quantity is also measurable as the luminosity distance d^ 



because of the reciprocity theorem of Etherington 10| (see also Ellis 14]), 



dL = {l + zyC{wo,y). (8) 

(iii) The Maximum of Observer Area Distance. Generally speaking, C{wo,y) reaches a 
maximum Cmax for a relatively small redshift Zmax (Hellaby 12]; see also Ellis and Tivon 



[13] and Araiijo and Stoeger [ij]). At Cmax, oi course, we have 

dC(wo,z) dC{wo,y) 



dz dy 

further conditioned by 

d^C{wo,y) 



0, (9) 



dz^ < '■ (1°) 



Furthermore, of course, as we shall review below, the data set will give us y = y{z), from 
which we shall be able to find ymax = ymax{zmax)- These Cmax and Zmax data provide 
additional independent information about the cosmology. Without Cmax and Zmax we cannot 
constrain the value of A. 

(iv) Galaxy Number Counts. The number of galaxies counted by a central observer out 
to a null radial distance y is given by 

rv 
N{y)=ATi ij{wo,y)m B{wo,y)C{wo,y) dy, (11) 

Jo 

where /z is the mass-energy density and m is the average galaxy mass. Then the total energy 

density can be written as 

dz 1 

fi{wo,y) = mn{wo,y) = Mo{z) — -— -, (12) 

dy B{wo,y) 

where n{wo, y) is the number density of sources at {wq, y), and where 

m 1 1 dN 
J dQ Tq dz 

Here dQ is the solid angle over which sources are counted, and J is the completeness of the 

galaxy count, that is, the fraction of sources in the volume that are counted is J. The effects 

of dark matter in biasing the galactic distribution may be incorporated via m and/or J . In 

particular, strong biasing is needed if the number counts have a fractal behaviour on local 



Mo^-^-,— (13) 



scales (Humphreys et al [l5j). In order to effectively use number counts to constrain our 
cosmology, we shall also need an adequate model of galaxy evolution. We shall not discuss 
this important issue in this paper. But, fundamentally, it would give us an expression for 
m = m{z) in equations flT^ and flT^ above. 

There are a number of other important quantities which we catalogue here for complete- 
ness and for later reference. 

First, there are the two fundamental four- vectors in the problem, the fluid four- velocity 
m" and the null vector /c", which points down the generators of past light cones. These are 
given in terms of the metric variables as 

^a ^ ^-l^a^ ^ },a ^ ^AB)-'5%. (14) 

Then, the rate of expansion of the dust fluid is 3H = VaW'', so that, from the metric (1) 
we have: 

8 



where a "dot" indicates d/dw and a "prime" indicates d/dy, which will be used later. For 
the central observer H is precisely the Hubble expansion rate. In the homogeneous (FLRW) 
case, H is constant at each instant of time t. But in the general inhomogeneous case, H 
varies with radial distance from y = on t = tg- From our central conditions above (3), we 
find that the central behavior of H is given by 

as.^O: i/(».,)^^|M=i/(».0). (16) 

At any given instant w = Wq along y = 0, this expression is just the Hubble constant 
Hq = H(wo, 0) = Aq Bq (-B)o as measured by the central observer. In the above we have 
also written Aq = A{wq, 0) and Bq = B{wo, 0). 

Finally, from the normalization condition for the fluid four- velocity, we can immediately 
see that it can be given (in covariant vector form) as the gradient of the proper time t along 
the matter world lines: u"" = —t,a. It is also given by ([1]) and flT^ as 



Ua = 9ahU^ = -Aw^a + By^a- (17) 

Comparing these two forms implies 

dt = Adw - Bdy ^ A = t^ , B = -ty, (18) 

which shows that the surfaces of simultaneity for the observer are given in observational 
coordinates by Adw = Bdy. The integrability condition of equation flTSl) is simply then 

A' + B = 0. (19) 

This turns out precisely to be the momentum conservation equation, which is the key 
equation in the system and essential to finding a solution. 

IV. THE SPHERICALLY SYMMETRIC FIELD EQUATIONS IN OBSERVA- 
TIONAL COORDINATES 

Using the fiuid-ray tetrad formulation of the Einstein's equations developed by Maartens 
|16l | and Stoeger et al [l7|, one obtains the Spherically Symmetric field equations in obser- 
vational coordinates with A 7^ (see Stoeger et al J5| for a detailed derivation). Besides the 
momentum conservation equation (IT9|) . they are as follows. 



A set of two very simple fluid-ray tetrad time-derivative equations can be quickly inte- 
grated to give: 

B{w,y) C^{w,y) 

where /i^ again is the relativistic mass-energy density of the dust, including dark matter, 
and uq, that is u specified on w = wq, is a quantity closely related to fimg (see equation (128|) 
below). In deriving and solving these equations, and those below, we have used the typical 
A equation of state, pa = — /^a, where pa and /ia = g;^ are the pressure and the energy 
density due to the cosmological constant. Both uq and /io are specified by data on our past 
light cone, as we shall show. fi\ will eventually be determined from the measurement of 
Cmax and Zmax- Essentially u is defined by the right-hand-side of equation (1211) . 
The fiuid-ray tetrad radial equations are: 



c — cyA ^ B J 2-° /^"i' y"^"^) 

(ujoiy) + if)C'{wo,y)\' = -i/i„.o B{wo,y) C\w,,y) ff + §); (23) 



% = %%-{u + ^)aB. (24) 

The remaining "independent" time-derivative equations given by the fiuid-ray tetrad formu- 
lation are: 



B B A , o 1 ,o / V 

rom equation fl23|) we see that there is a naturally defined "potential" (see Stoeger et al 
5|]) depending only on the radial null coordinate y - since the left-hand-side depends only 
on y, the right- hand- side can only depend on y: 

10 



where N^,{y) is an arbitrary function, whose central behavior is the same as that of number 
counts (Stoeger et al |I5{]). Thus, 

uJo{y) = -^ - 2C\w, y) /^-°(^) ^^^o,y) C\w,,y) F{y) dy. (28) 

Connected with this relationship is equation (|2T1) . which we rewrite as 
CC A C'^ AB AB 







(29) 



CC 2B C^ 2^2 C^ 

where Cq = C{wo, y). 

Stoeger et al ^] and Maartens et al [IS"] have shown that equations flTTI) and fl2^ can be 
transformed into equations for A and 5, thus reducing the problem to determining C: 

C 

^ ^ [F2 - 1 - 2(.;o + ^^K/Q)Cl/C + (/iA/3)C2]i/2 ^^^^ 

C" 
^ ^ F ± [F2 - 1 - 2(^0 + /iA/6)CoVC + (/iA/3)C2]i/2 ■ (^^) 



The Lemaitre-Tolman-Bondi (LTB) exact solution (Lemaitre [19]), Tolman 20|, Bondi 21| ; 
and cf. Humphreys |22] and references therein) is obtained by integration of (!30!) along the 
matter flow y =constant using flTSl) 

^ " ^^y^ = y [F2 - 1 - 2(a;o + I1k/Q)CI/C + (/iA/3)C72] V2 ' (^2) 

where T(|/) is arbitrary, provided we identify 

F^ = l-kf, k = 0,±l. (33) 

Here / = /(?/) is a function commonly used in describing LBT models in the 3 + 1 coordi- 
nates (Bonnor 23|). 



V. INTEGRATION WITH FLRW {k = AND A / 0) DATA 

In this section and the next we use a generalization (to incorporate the Cmax and z„ 
data) of the integration procedure described in detail in Araiijo ei a/ [2|, which in turn is 



•^max 







an improvement of the integration scheme developed by Araiijo and Stoeger p|, to solve the 
above system of SS field equations when A 7^ 0. First, we consider a concrete, simplified, 
but very relevant example. Suppose that we know that the universe is flat. Then F{y) = 1. 
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This means that we need only the observer area distance, or the galaxy number counts - 
not both. F = 1 establishes a relation between these data functions. Suppose also that, 
though we do not know that the universe is FLRW, we find that our observer-area-distance 
and galaxy-number count data can be fit - or very closely approximated - by the FLRW, 
A 7^ observational relationships as functions of the redshift z. For the fiat case these are 
(Araiijo and Stoeger |l4|): 



'nAVSHoiil + V^)Qiz) - (1 - V3)] I 



Q{z),K 



Cn 



-1 



QiO),K 



(l + v^) + (l + z) ,3/3131' 



A 



SHo 



21 



/2 + y3 

K = \ ■ , 



(34) 



and 



AoHoJ[nA + {l + zf{l-QA)] 



Cn~^{u, n) is the inverse of the Jacobi elliptic function Cn{u,n), where n is the modulus. 
Equation (!34l) can clearly also be written in terms of elliptic integrals of the first kind 
(Araiijo and Stoeger |l4|). 



Equations flMl) and fl35|) are the A 7^ analogues of the familiar characteristic FLRW 



Tj) = rQ{z) and Mq = Mq{z) relationships for A = (Ellis and Stoeger |2j]; Stoeger, et al. 
[51]): If the universe is FLRW and A = 0, tq^z) and Mo(z) will have those functional forms. 
Equations flM|) and fl35|) - or their elliptic-integral equivalents - are the corresponding 
characteristic functional forms for fiat FLRW, A 7^ cases. As in the A = cases, however, 
it is not a trivial conclusion that data which satisfies equations flM|) and fl35l) implies an 
FLRW, A 7^ universe. This must be demonstrated - and can be demonstrated - by using 
equations (1531) and fl5S]) as data functions to solve the field equations and obtain an FLRW 
solution. This was done for the fiat A = case by Araiijo and Stoeger [l|. We now do the 
same for these fiat A 7^ cases. 
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Because of the additional parameter (A) in the equations, though, in the general (non- 
flat) case we need more data than supplied simply by equations flM|) and fl35l) . When we 
know the universe is flat, this extra data can be supplied by equation (1351) . which will then 
automatically enable us to calculate ^a - from the flatness condition. But also on w = Wq 
we can observationally determine where Cq = tq reaches its maximum value Comax and the 
redshift Zmax at which this occurs. These measurements provide the needed extra data in 
the general case, and can also be used, instead of the data in equation fl35|) in the flat case. 
From equation (l34l) . we can immediately determine the equation for this maximum redshift, 
which will be 

dro{z)/dz = 0. (36) 

Plugging the observationally determined values of Zmax into this equation, we obtain a 



unique relationship between Zmax and ^a (Araiijo and Stoeger [IJ]), since Hq cancels out of 
equation (136|) . Using this relationship along with Comax in equation (1M|) will also determine 
Hq. The precise interpretation and deflnition of these parameters, e.g. that ^a ^ depending 
not only on A but also on Hq, the FLRW Hubble parameter - is the density parameter 
for A, is in reference to a supposed FLRW universe {Hq, and therefore ^a, cannot be 
unambiguously deflned in a general exactly spherically symetric - also often referred to as 
a Lemaitre-Tolman-Bondi (LTB) - universe; the deflnition of the rate of expansion given in 
equation flTSj) is not the only one that could be chosen, or measured). That assumption is 
validated by continuing with the integration and showing that a universe with such data is 
indeed FLRW. 

Of course, this is a somewhat contrived case, since we would not attempt to flt the data 
to such a functional form flM|) unless we already suspected that the universe may be FLRW, 
and that therefore Qa and Hq, the FLRW density parameter and the Hubble parameter 
at w = Wq and y = 0, can be deflned in terms of an FLRW model. But besides being 
a very relevant simple case, it serves to illustrate the integration procedure with deflnite 
meaningful data input functions. 



0.' 



Solving the null Raychaudhuri equation (!22l) on w = Wq with this data (see Stoeger et al 
and also Araiijo and Stoeger [sl]) yields the following relation between redshift and the 
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null coordinate y: 

(1 - V3) - (1 + ^/?,)Cn{Ly + a) 



l + z 



■'1^. - 1 



Cn{Ly + cr) - 1 



L = -^Qj,V3AoHo{ — - 1 = , 

V i'A UJO 



a = Cn~^ 



g(0),«: 



(37) 



Using equation ([6]) we can now write A[wq^ y) as a function of y as 



^^'^^^--Mk- %-J^:^cL..^ ,38) 

Aq = A{wo, 0) is a constant scaling factor which can be chosen arbitrarily. 



If we wish to continue using our assumption that this will be an FLRW universe, then 
the observer area distance C{wq, y) as a function of y can be written as 



If we don't wish to assume FLRW here, we can still find C{wq, y) by using equation (l3l|l in 
conjunction the result given in equation ( l37j) . 



Furthermore, we can clearly now determine what ymax is, corresponding to z^axi from 
equation (1371) . 



Now we begin to move our solution off our past light cone, w = wq. Since y is chosen to 
be a comoving radial coordinate, the functional dependence of A{w, y) with respect to y can 
not change as we move off our light cone. We have already mentioned the freedom we have, 
temporally setting aside central-condition considerations, to rescale the time coordinate 
w, which is affected by choosing A{w,Q). Therefore, this freedom effectively corresponds 
to choosing the functional dependence of A{w,y) with respect to w in any way we like, 
constrained only by the form of A{wo, y) [later this choice may have to be adjusted to satisfy 
all the central conditions, those on C{w,y) and B{w,y)]. In our expression for A{wo,y) is 
hidden an implicit dependence on w. We need to extract that dependence and make it 
explicit, so that we can then determine the general dependence of A on w and proceed with 
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the integration. In general, this is not simply achieved by replacing wq with w, because - 
besides the wo dependence arising from setting w = wq when we write equation ([6]) - there 
is another part of the wq dependence which derives from integration constants of the null 
Raychaudhuri equation and remains through the entire problem. 

At this point, accordingly to what we just pointed out, we arbitrarily set the w dependence 
for A and proceed with the integration. The next step is then the solution of equations (fT9|) 
and fl30|) to determine B and C respectively. Their formal general solutions are: 



B = - A'dw + l{y), 
where l{y) is determined from the condition A{wo,y) = B{wo,y), and 



(40) 



C 



^-^ + l]CoMnh^'l^^' 



Adw + h{y) 



1/3 



(41) 



fJ'A J VV 4 

where h{y) is determined from the data tq = C{wq, y). Since we know C{wo, y) from equation 
( l39|) . or from equation ( l34l) with equation (1371) . u!o{y) is obtained from equation (1311) and is 
given by 



UJo 



2^2 



1 



:i+z)c' 



n2 



where we have used (13]) and (El) to write 



B{wo,y) = A{wo,0)/[l + z 



(42) 



(43) 



Here z{y) is given by equation (1371) . When we do not know that the universe is flat, we 
must, of course, first determine C{wQ,y), in order to determine F{y) from equation (127 ). 
This can be easily done, as explained in Maartens, et al |l8|, in Araiijo and Stoeger |l4| . 
and in the next section. 



B{w,y) and C{w,y) are then determined by integrating equations (lT9l) and (l30l) with 
respect to w. B{w,y) and C{w,y) are further constrained, as discussed above, by the fact 
that they have to satisfy the central conditions ([5]). Now, it is clear from an examination 
of these equations that unless A{w,y) has a very specific functional dependence on w the 
resulting solutions B{w,y) and C{w,y) will not satisfy the central conditions. That implies 
that, although we can find solutions to the field equations, it does not guarantee that the 
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null surface on which we assume we have the data is a past lisrht cone of our world line (Ellis 

fl 
et al |4|). So we conclude that given the fulfilment of the following conditions: 

(1) A{wQ^y) is determined by the data and the central conditions; 

(2) The coordinate y is choosen to he a comoving radial coordinate; 

(3) The central conditions (^; 

we can remove the freedom of rescaling the time coordinate w and completely determine 
A{w, y). Thus, all the coordinate freedom in y and w has been used up at this stage. 
Therefore, following the above analysis we find that the appropriate form for A{w,y) is 



A{w,y) = -Ao 



1 



1 



1-Cn 



fa - 1) + ^^ 



\ V ^A 7 (1 - V3) - (1 + V3)Cn 
where z{y) is given by equation fl37|) . Now, observing that 

dz 
dy 



V-l) + ^^ 



AoHoJnj, + {l + zr{l-n,,), 



(44) 



(45) 



we substitute A'{w,y) and A{w,y) into equations (HOj) and (HTj) and determine the arbitrary 
functions of y that arise from these integrations by the conditions B{wQ,y) = A{wQ,y) and 
C{wo,y) = Tod/) respectively. Thus, 



B{w,y) = A{w,y) 



-An 



n^ 



and 



1 



1-Cn 



fa - 1) + ^^ 



(1 - v^) - (1 + V3)C 



n 



a 



C{w, y) = A{w, y)y = -Aoy ( fl ^ " 1 



1-Cn 



(a-l) 



too 



w-y 

WQ 



(46) 



:^-i) + 



w-y 



(47) 



;i - V3) - (1 + V3)Cn 

which are the FLRW form of the solutions for A 7^ in observational coordinates. One can 
easily check (after some algebra) that the central conditions ([5]) are all satisfied, which in 
turn guarantees that the null surface on which we assume we have the data is indeed a past 
light cone of our world line. 



VI. THE GENERAL SOLUTION IN THE NON-FLAT CASE 

We now outline the integration procedure in the case where we do not know whether 
the universe is fiat or not, and where data gives us redshifts z, observer area distances 
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(angular-diameter distances) tq^z), "mass source densities" Mo{z) which cannot be fit by 
the FLRW functional form, and the angular-distance maximum Cmax{wo,z) at Zmax- It is 
important to specify the latter, because, as we have already emphasized, without them, 
we do not have enough information to determine all the parameters of the space-time in 
the A 7^ case. For instance, although we can determine C{wq, z) with good precision (by 
obtaining luminosity distances di and employing the reciprocity theorem, equation ([8])) 
out to relatively high redshifts, at present we do not yet have reliable data deep enough 
to determine Cmax and z^ax- But this has just recently become possible with precise 
space-telescope distance measurements of distance for supernovae la. 

In pursuing the general integration with these data, we use the framework and the 
intermediate results we have presented in Section IIVI Obviously, one of the key steps we 
must take now is the determination of the "potential" -F(|/), given by equation (1271) . This 
was done in a similar way for A = by Araiijo and Stoeger [l|, as indicated above. This 
means we need to determine C'{wo,y) and C{wo,y), which we now write as Cq and Cq, 
respectively. We also need A{wo,y). We remember, too, that at on w = wq we have chosen 
B{wo,y) = A{wo,y), which we have the freedom to do. 

Clearly, Cq can be determined from the r'o(z) = C{wq,z) data, through fitting, along 
with the solution of the null Raychaudhuri equation fl22|) . as indicated in Section |Vl to 
obtain z = z{y). A{wo,y), too, is obtained from redshift data along with this same z{y) 
result. Cq is somewhat more difficult to determine. But the procedure is straight-forward. 

We determine Cq by solving equation (^^ for it on w = wq. Using equations ([3]) and 
( TT9|) . we can write this now as: 



C^iy) _ A'Qiy)C'Qiy) ^, 



Co{y) Ao{y)Coiy) 



Af^{y){uJo{y)+i2j2). (48) 



But, from equation (!2T!) we can write u!o{y) in terms of Co{y), C'^iy), and C*o(l/)- So equation 
148) becomes: 

• C^(y)CUyl Al{y) A'^iy) {C',{y)f , Al{y)Co{y) 
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This is a linear differential equation for Co{y), where from data we know everything 
on our past light cone, w = Wq, (once the null Raychaudhuri equation fl22l) has been 
solved) except Co{y) itself and /ia, which is a constant that can be carried along and 
determined subsequently from C{wo,Zmax) and z^ax measurements (see below). Thus, we 
can easily solve equation fH9|) for Coiv), which will also depend on the unknown constant fi\. 



However, introducing this result back into equation fH9|) . and evaluating it at Umax, which 
corresponds to Zmax, we have simply 

'^"^"-^ = ^oh::^) 2 *"- e^"' 

where Coiymax), as we have already emphasized, also depends on the unknown /ia- Since 
everything else is now known, equation fl50|) is now an algebraic equation for fiAiVmax), or 
equivalently for jiK^Zmax)- Obviously, we could have simply worked out this result in terms 
of Zmax to begin with. 

With this determination of /^a, we know CQ{y) completely, and can now determine F{y) 
from equation (1271) . From there on, we can follow the solution oS w = wq for all w, as we 
have outlined in Section |Vl Obviously, we shall obtain very different results than we did 
there for fiat FLRW data - depending on the exact character of our more general data. 
This completes the framework for solving these exact spherically symmetric field equations 
for adequate data when A 7^ and the space-time is not fiat. 



VII. CONCLUSION 

In this paper we have shown in detail how to construct fiat dust-filled A 7^ Friedmann- 
Lemaitre-Robertson- Walker cosmological models from FLRW cosmological data on our 
past light cone, by integrating the exact spherically symmetric Einstein field equations in 
observational coordinates, extending the approach developed by Araiijo and Stoeger pj 
and improved in Araiijo et al 2|. Besides being of interest in its own right - demonstrating 
how such data fully determines the models - it also illustrates in a simple case how the 
more general SS Einstein equations can be integrated in observational coordinates with 
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data arrayed on our past light cone, w = Wq- Then, we have gone on to show how to 
integrate these exact spherically symmetric (LBT) equations, also for A 7^ in cases where 
the data are not FLRW, and the space-time is not known to be flat. It is essential for 
these to have data giving the maximum of the observer area (angular-diameter) distance, 
Co{wq, -Zmax)? and the redshift z^ax at which that occurs. This enables the determination of 
the vacuum-energy density fi\, which would otherwise remain undetermined. 
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